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We have studied the origin of a counter intuitive diffusion behavior of Fe and N atoms in a iron 
mononitride (FeN) thin film. It was observed that in-spite of a larger atomic size, Fe tend to diffuse 
more rapidly than smaller N atoms. This only happens in the N-rich region of Fe-N phase diagram, 
in the N-poor regions, N diffusion coefficient is orders of magnitude larger than Fe. Detailed self- 
diffusion measurements performed in FeN thin films reveal that the diffusion mechanism of Fe and 
N is different - Fe atoms diffuse through a complex process, which in addition to a volume diffusion, 
pre-dominantly controlled by a fast grain boundary diffusion. On the other hand N atoms diffuse 
through a classical volume-type diffusion process. Observed results have been explained in terms of 
stronger Fe-N (than Fe-Fe) bonds generally predicted theoretically for mononitride compositions of 
transition metals. 


I. INTRODUCTION 

Transition metal nitrides (TMN), specially 3d TM 
mononitrides (TMMN) are an important class of mate¬ 
rials exhibiting several interesting properties such as su¬ 
per hardness^”^, superconductivity^, corrosion and wear 
resistance^’^, etc. Along the 3d series, there is a char¬ 
acteristic variation in the heat of formation (Ai7p for 
the 3d TMMN - an initial increases to a maximum is 
followed by a decrease and a plateau in AH^.^ This in¬ 
herently makes formation of early 3d TMMN viz. ScN, 
TiN, VN, GrN easier than that of late ones viz. MnN, 
FeN, GoN, NiN; e.g. Ai7j=-338kJmol“^ for TiN and 
-47kJmol“^ for FeN at 298 K. As a result, mononitrides 
of Ti,V,Cr can be easily prepared and possess excel¬ 
lent thermal stability due to which they have been in¬ 
tensely investigated.On the other hand, magnetic 
mononitrides (e.g. MnN,FeN, GoN) started to gain at¬ 
tention rather recently.Unlike early 3d TMMN, FeN 
or GoN can only be formed in the form of thin films 
using non-equilibrium processes such as reactive sputter- 
i]^gi4,i5,22,23^ pulsed laser deposition^"^, and more recently 
by molecular beam epitaxy assisted with a rf-discharge 
nitrogen/ammonia source.Nickel mononitrides are 
yet to be evidenced experimentally.On the basis of 
energetics of mononitrides, FeN and CoN are expected 
to be metastable. 

The metastable nature of FeN turned out to be a boon 
as FeN films were exploited as a source of spin-injection 
to semiconductors or diluted magnetic semiconductor in 
spintronics.^^ FeN when heated above 650 K yields a ther¬ 
mally stable 7'— Fe 4 N giving rise to an array of litho¬ 
graphically defined spin-valves.The mechanism lead¬ 
ing to such structural transformation was assumed to be 
controlled by N-diffusion. XPS measurements for Fe 2p 
and N Is peaks were used to measure N diffusion in Fe.^^ 
However, the conclusion drawn from such measurements 
that fast N diffusion leads to such transformations can 


be misleading since interdiffusion of N in Fe, not self¬ 
diffusion of N was measured. 

Recent Fe and N self-diffusion measurements per¬ 
formed using neutron reflectivity (NR) show that N self¬ 
diffusion is slower than Fe.^^ This is a counter-intuitive 
result, defying established diffusion models for binary 
metal-metalloid systems where a smaller atom always 
diffuses faster than a larger atom.^^ In absence of a suit¬ 
able radioactive tracer, N self-diffusion is rather difficult 
to measure.Nevertheless by using ^^N labelling, it 
can be obtained by doing depth-profile measurements 
using secondary ion mass spectroscopy (SIMS)^^’^^, nu¬ 
clear reaction analysis (NRA)^ and NR. Among these, 
later provides an unique opportunity to measure self¬ 
diffusion lengths down to 0.1 nm and a possibility to mea¬ 
sure N self-diffusion in the low temperature regime (be¬ 
low 500 K).i^’34 

In the present work, we have carried out a study of 
Fe and N self-diffusion process to understand the origin 
of anomalous diffusion in non-magnetic FeN compound. 
Thin film samples were prepared using nitrogen alone as 
a sputtering gas in a reactive magnetron sputtering pro¬ 
cess. X-ray diffraction, absorption and Mossbauer spec¬ 
troscopy measurements confirmed ZnS-type structure of 
samples. NR measurements show that N diffusion is 
slower than Fe, however detailed diffusion mechanism was 
obtained from SIMS depth profiles using Le Glaire’s anal- 
ysis35,36 grain-boundary {gb) diffusion. It was found 
that in the low temperatures regime (up to 550 K) fast 
Fe diffusion takes place predominantly through the gh 
regions while N diffusion is a conventional volume type 
diffusion. As the temperature is increased beyond it, the 
difference between Fe and N diffusion decreases leading 
to structural transformations - essentially triggered by Fe 
diffusion and followed by N diffusion to the extent that 
N diffuses out of the system. This is an important result 
defying a general misconception the fast N diffusion leads 
to structural instability and can be applied to understand 
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the thermal stability of transition metal mononitrides. 


II. EXPERIMENTAL 

In a direct current-magnetron sputtering (dc-MS) tech¬ 
nique an iron target (purity 99.95%) was sputtered us¬ 
ing nitrogen (purity 99.9995%) alone as the sputtering 
medium at a constant power of 100 W. N 2 gas was flown 
at a constant flow of 10 seem yielding a pressure of about 
0.4 Pa, while the pressure before gas flow was about 
lx 10“^ Pa. Fe target was pre-sputtered for 10minutes 
using Ar gas to remove surface contaminations. Follow¬ 
ing samples were prepared at ambient temperature on 
Si(lOO) and float glass substrates: 

(Nl): [FeN(7.5 nm) |^^FeN(7.5 nm)] xio 
(N2): [FeN(7.5nm)|Fe^^N(7.5nm)]xio 
(S): [FeN(100 nm)|^^Fe^^N(2 nm)|FeN(100 nm)] 

Samples (Nl) and (N2) were used to measure Fe and N 
self-diffusion using NR, sample (S) was used to measure 
Fe and N self-diffusion simultaneously using SIMS. Nat¬ 
ural Fe and ^^Fe (^95% enriched) targets were sputtered 
using natural nitrogen and ^^N (^98% enriched) gases. 
To avoid any mixing of nitrogen isotope gases, chamber 
was evacuated after deposition of each layer and gas flows 
were monitored using a residual gas analyzer (RGA). 

The long-range structure transformation in the sam¬ 
ples were studied using x-ray diffraction (XRD) using a 
standard diffractometer (Bruker DSAdvance) equipped 
with CuK-a x-ray source and a silicon stripe detector 
(Lynxeye). The local structure was investigated using x- 
ray absorption spectroscopy (XAS) and conversion elec¬ 
tron Mossbauer spectroscopy (GEMS). The XAS mea¬ 
surements at Fe K-edge were performed at BL-9 beam- 
line^^ and those at N K-edge at BL-1 beamline^^, both 
at Indus-2 synchrotron radiation source at Indore. While 
Fe K-edge measurements were carried out in fluorescence 
mode, those at N K-edge were measured in total electron 
yield mode. To study the phase transformation, sam¬ 
ples were annealed using a vacuum furnace for about two 
hours at each temperatures. NR measurements were per¬ 
formed at AMOR reflectometer at SINQ/PSI, Villigen, 
Switzerland. SIMS measurements were carried out on a 
Hiden Analytical SIMS Workstation using O J as primary 
with 5keV energy and 400 nA beam current. The base 
pressure in SIMS measurement chamber was lx 10“^Pa 
and during measurement pressure was of the order of 
lx 10“^Pa due to differentially pumping of oxygen gas 
source. 


III. RESULTS AND DISCUSSION 

A. Structural measurements 

Figure 1(a) shows XRD patterns of sample (Nl) in the 
as-deposited state and after annealing at various tem¬ 
peratures. In the as-deposited state, a single broad peak 



FIG. 1. XRD (a) and CEMS(b) patterns of sample (Nl) in 
the as-deposited state and after annealing at different tem¬ 
peratures. 


appears around 20=56 degree. Erom GEMS and XAS 
measurements (shown later), it was confirmed that the 
formed phase is y'^-PeN having ZnS-type structure ori¬ 
ented along (220) plane.After annealing at 475 K, a 
shift in the peak position towards higher angle was ob¬ 
served. However, the XRD pattern of sample annealed at 
525 K is almost identical to that of previously annealed 
sample. Such shift in the peak position is an indication of 
annihilation of free volume leading to densification of the 
structure.It is interesting to note that the average 
crystallite size (qs) (calculated using Scherrer formula^^) 
remains at a value of about ^ 5 nm upto an annealing 
temperature of 525 K. At a higher temperature of 575 K, 
additional peaks corresponding to e — Ee 3 _zN( 0 <z<l) 
start to appear. With further increasing annealing tem¬ 
peratures, growth of 5 —Ees-zN phase can be seen. Along 
with it a continuous shift of peaks position towards higher 
20 values were observed indicating % ^ O’. 

Eigure 1 (b) shows GEMS spectra for sample (Nl) in 
the as-deposited state and after annealing. Obtained 
GEMS spectra were fitted using a computer program 
NORMOS/SITE^^ and fitted parameters are given in ta¬ 
ble I. As-deposited and the sample annealed at 525 K 
shows an asymmetric doublet which is typically observed 
for 7 '"-EeN in ZnS-type structure. Such spectra can be 
de-convoluted into two singlets: one with a smaller and 
other with a larger value of isomer shift (S). The singlet 
with (5=0.01 mm/s corresponds to Ee surrounded tetra- 
hedrally to N, while other singlet originates due to de¬ 
fects and vacancies.Annealing above 575 K results 
in appreciable changes in GEMS spectrum, at this tem¬ 
perature it can be best fitted using mixture of two sin¬ 
glets and two quadrupole split doublets corresponding to 
7 '" and 5 — Ees-zN phase, respectively. Obtained fitting 
parameters match-well with the reported values.It 
is known that 5 — Ees-zN phase exist in a wide composi- 
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TABLE L Conversion electron Mossbauer spectroscopy 
(GEMS) parameters for iron mononitride thin films in the 
as-deposited state (300 K) and after annealing at various tem¬ 
peratures. Here S stands for a singlet, Q for a doublet and H 
for a sextet, S for isomer shift, E for quadrupole splitting, H 
for hyper fine field and RA for relative area. 


T 

Component 

(5 

r 

H 

RA 

(K) 

( 

mm/s) 

(mm/s) 

(Tesla) 

(%) 



±0.03 

±0.03 

±0.2 

±2 

300 

SI 

0.10 

- 

- 

54 


S2 

0.57 

- 

- 

46 

525 

SI 

0.11 

- 

- 

57 


S2 

0.60 

- 

- 

43 

575 

SI 

0.11 

- 

- 

40 


S2 

0.60 

- 

- 

34 


Q1 

0.46 

0.31 

- 

23 


Q2 

0.30 

0.50 

- 

3 

625 

HI 

0.38 

0.0 

8.7 

13 


H2 

0.31 

0.0 

17.4 

17 


H3 

0.33 

0.65 

- 

70 

675 

HI 

0.40 

0.0 

7.6 

22 


H2 

0.31 

0.0 

21.7 

68 


H3 

0.30 

0.0 

3.6 
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tion range in which its magnetic properties also get tuned 
with nitrogen concentration.It was observed that, at 
room temperature, as 'z e — Fes-zN phase becomes 
non-magnetic, whereas as ‘ 2 ; ^0’ it becomes ferromag¬ 
netic. Appearance of sextet confirms ferromagnetic or¬ 
dering which is also supported by XRD results indicating 
‘z ^0’ in s — Fcs-zN. Obtained fitting parameters (ta¬ 
ble I) for GEMS spectra measured above 575 K match- 
well with the reported values. 

To get precise information about the local structure of 
Fe and N, XAS measurements at Fe and N K-edge were 
performed on sample (Nl) and are shown in figure 2. 
A strong pre-edge peak around 7113 eV can be seen be¬ 
fore Fe K-edge. This is a signature of quadrupole transi¬ 
tion and its probability strongly depends on the direction 
of the local electric field, which in turns gets influenced 
by local site symmetry.For such a pre-edge peak to 
appear, inversion symmetry must be broken (transition 
from a bound core level to a higher level empty state). 
This is only possible when Fe is surrounded tetrahedrally 
to N atoms because this arrangement is asymmetric un¬ 
der inversion.This clearly indicates ZnS-type structure 
(tetrahedral coordination of Fe atoms) in our samples. 
Such pre-edge feature has been frequently used to as¬ 
sign tetrahedral or octahedral coordination in transition 
metal complexes such as Mn^^, Cr^^, etc. To further con¬ 
firm this result, XAS spectra was also taken at N K-edge 
(shown in the inset of figure 2). Apart from edge feature 
‘a’ at 401 eV, three features that are assigned as (c, d and 
e), are observed due to the electronic transition to empty 
N-2p state hybridized to Fe-3d state. Various TMMN 
(viz. Ti, Cr, V) having NaCl type structure, the feature 
‘a’ splits into two components due to crystal field split- 



Energy (eV) 

FIG. 2. Fe K-edge XAS spectra of iron mononitride thin films. 
Inset of the figure shows spectra taken at N K-edge. 

ting caused by the octahedral coordination of N atoms 
surrounding the metal ions.^^’^^ Absence of any splitting 
around the feature ‘a’ is a clear indication of tetrahe¬ 
dral coordination of Fe surrounding N atoms. Less in¬ 
tense features (b) and (d') observed in the spectra are due 
to surface oxidation state. Combining the information 
obtained from XAS measurements, it can be confirmed 
that FeN phase has a ZnS-type structure. This is also 
in agreement with low temperature high field Mossbauer 
spectroscopy measurements in this compound. 

B. Self-diffusion measurements 

NR is a precise tool to measure atomic self-diffusion 
and diffusivity as low as lxl0“^^m^s“^ have been mea¬ 
sured using this technique.It is known that 
neutron scattering length {bn) varies for isotopes for 
naturalp^^ natural^^ 15^ 6n=9.45fm, 2.3 fm, 

9.36 fm and 6.6 fm, respectively. Therefore, periodic iso¬ 
tope multilayer are widely used to study atomic self¬ 
diffusion using Figure 3 (a) and (b) 

shows NR patterns for samples (Nl) and (N2), respec¬ 
tively. Bragg peaks appearing due to isotopic contrast of 
natpg/57pg and ^^^N/^^N can be seen clearly. Patterns 
were fitted using a computer program based on Parratt 
formulism^^’^^ and obtained layer thickness for sample 
(Nl) is 7.9 nm, (N2) is 7.3 nm, close to nominal values. 

After annealing the intensity of Bragg peak start to 
decay and such decay is more rapid for sample (Nl) than 
for (N2). A decay of Bragg peak intensity is a measure 
of atomic diffusion across interfaces. Obtained results 
clearly indicate that Fe self-diffusion is faster than N. 
Detailed fitting of NR data yields diffusion length(Ld), 
which are plotted in the inset of figure 3(b). Clearly, Ld 
is significantly larger for Fe than for N. This result, al¬ 
though counter intuitive, agrees well with previous stud¬ 
ies on 7 '"-FeN compound.To get further insight leading 
to such anomalous behavior, we did detailed SIMS mea¬ 
surements on sample (S). With SIMS, unlike NR, both 
Fe and N self-diffusion can be measured simultaneously. 

Figure 4 shows SIMS depth-profile in sample (S) for 
^^Fe(a) and ^^N(b), isochronally annealed for 1 hr at each 
temperature. As annealing temperature is increased, the 
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FIG. 3. NR patterns of sample (Nl) (a) and (N2)(b) in the 
as-deposited state (ASD) and after annealing at different tem¬ 
peratures for 2 hours. Here scatter points are experimental 
data and solid line is ht to them. Inset of hgure(b) shows 
variation of diffusion length (Ld) obtained by htting. Here 
solid and open symbols represent for Fe and N diffusion, re¬ 
spectively. Typical error bars in calculating Ld are less than 
size of symbols. 
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FIG. 4. SIMS depth prohle of ^^Fe(a) and ^^N(b) for sample 
(S) after annealing at different temperatures. 
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FIG. 5. SIMS depth-prohle for sample (S) showing ^^Fe(a) 
and ^^N(b) distribution plotted versus to calculate volume 
diffusion, in the as-deposited and annealed state. Inset of (a) 
shows depth-prohle of ^^Fe plotted versus to calculate 
grain-boundary diffusion. 


width of ^^Fe profiles become broader than those of 
indicating faster Fe self-diffusion than N, in agreement 
with NR results. For such SIMS depth-profiles, concen¬ 
tration profile of an element at a depth z can be fitted 
using a Gaussian distribution function to obtain volume 
diffusivity(I)y) using^^: 




const. ( —z^ \ 

\4D7t) 


( 1 ) 


Here t is annealing time at a temperature T. Alterna¬ 
tively, a linear relation between ln{c) and z^, yields slope 
of the fitted line, given by: slope= —l/4Dyt. Obtained 
N diffusivities are shown in 6(c). While ^^N profiles fits 
well to this equation, ^^Fe profiles cannot be fitted using 
equation 1 alone, specially above 475 K. As shown more 
clearly in figure 5, ^^N profiles have single slope, ^^Fe pro¬ 
files have two distinct slopes. Fe diffusivity upto 475 K 
can be obtained using equation 1, above it, the additional 
slope can be fitted using Le Claire’s analysis^^ for gh dif¬ 
fusion. An skewness in ^^Fe depth-profiles above 475 K 
(marked by an arrow in figure 4(a)) is an indication of 
competing diffusion processes taking place through grain- 
boundaries in addition to a volume type diffusion. 
Since Dgi)»Dy^ it can be assumed that fast Fe diffu¬ 
sion observed in our samples is primarily gh diffusion. 
In this case, Dgij can be obtained from the slope of 4nc’ 
versus curve, using Suzuoka’s instantaneous-source 

solution equation^^’^^: 


sSDab = 1.084 


.^y/1.94 -5/2.91 


Here, s is segregation factor, S is width of gh. As such 
experimental estimation of 6 is difficult, a good approxi¬ 
mation is S ^0.5 nm, as often found in literature.Eq. 2 
is applicable for the condition in which a dimensionless 
quantity /3 = , has values between 10 and 100. 

Since in our case we find /3 ^ 23 at 500 K, equation 2 
can be applied to calculate sDgb. Inset of the figure 5(a) 
shows a linear fit to a curve plotted between 4nc’ and 
z^/^’ after annealing at 525 K. Similar analysis was car¬ 
ried out at other annealing temperatures and obtained 
values of sDgb and Dy are shown figure 6 (a) and (b), 
respectively. 

Grain-boundary diffusion is a complex process in which 
several competing processes may take place simultane¬ 
ously.^^ However for a given temperature range and du¬ 
ration, only few of them may get activated. Thus differ¬ 
ent kinetic regimes of Dgb have been observed, ascribed 
by Harrison.According to this classification, Dgb has 
three regimes known as type-A, B, and C. These regimes 
can be differentiated according to a relation between gh 
penetration depth {Lgb = (Dg^t)), volume penetration 

depth {Ly = (Dyt)) and grain size (gs)-^^ For type-C 
regime Lgb < gs, for type-B regime Ly < gs, and for 
type-A regime Ly > gs^^• Using these inequalities, ob¬ 
tained values Dy and Dgb can be divided into distinct 





















































5 


425 450 475 500 525 550 575 600 
^ 6x10 '^ 

CN * 

^ 4x10""’ 

Q 2x10'" 


9x10’^“ 

6x10'^° 

CO 

3x10'^° 

H—> 

;> 

^ 9x10'"° 
^ 6x10 "° 
3x10"° 


425 450 475 500 525 550 575 600 
Temperature (K) 

FIG. 6. Variation of Dg^ of Fe(a), Dy of Fe(b) and Dy of 
N(c) at different annealing temperatures. Typical error bars 
in estimation of diffusivity are of the order of size of symbols. 
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FIG. 7. SIMS depth-profile of ^^Fe(a) and ^^N(b) for sample 
(S) annealed at 525 K for different annealing times. Obtained 
values of Dy for Fe and N annealed for different times at 
525 K (c). Inset of figure (c) shows variation in Dgh of Fe. 
Arrhenius behavior of volume and grain-boundary diffusion 
of Fe and N (d). 


kinetic regimes as shown in figure 6. We find that for Fe 
diffusion, Lgi, < Qs {Lgi) ^2nm, gs ^5nm) for T<475K; 
for 475K<T<550K, Ly - gs and Ly > gs for T>550K 
corresponding to type-C, B and A kinetics, respectively. 

On the other hand for N diffusion, we find Ly < gs 
below 550K and above it Ly > gs. This indicates N 
diffusion is taking place within a grain below 550 K and 
as temperature is raised beyond it, multi-grain N dif¬ 
fusion takes over. It appears that Fe and N diffusion 
mechanism is significantly different. While only Dy takes 
place for N, Fe diffusion process seems to be more com¬ 
plex having distinct kinetic regimes C, B and A as an¬ 
nealing temperature increases. It is known that in type- 
C regime, diffusion takes place pre-dominantly through 
grain-boundaries leading to segregating of Fe atoms in 
the gh region. In type-B regime, in addition to gh diffu¬ 
sion, Dy starts. Though it is not too effective to cause 
any structural or magnetic transformation in this case. 
Finally above 550 K, Dy of both Fe and N becomes com¬ 
parable. Observed diffusion behavior gives a direct corre¬ 
lation between Fe and N self-diffusion with phase trans¬ 
formation. This is contrary to general preassumption 
that only N diffusion leads to phase transformation. 

Since isochronal diffusion measurements only give a 
snap shot of diffusion process, more insight about in¬ 
volved diffusion mechanism was obtained by doing de¬ 
tailed isothermal diffusion measurements between 475- 
550 K (for Fe) and 525-600 K (for N) in the steps of 
25 K for different annealing times. Representative SIMS 
depth-profiles taken at 525 K after various annealing 
times are shown in figure 7 for Fe (a) and N (b). Fol¬ 


lowing a similar process Dy and was calculated and 
shown in figure 7(c). Time-dependent diffusivity data 
can be fitted using: 

D = Dr + A • exp{—t/T) (3) 

Where, Dr is diffusivity in relaxed state, A a constant, 
t annealing time, and r is relaxation time. Using 3 we 
get r =1445(±300)s, and 2192(±300)s, respectively for 
Fe and N Dy; while for Fe r=2986(±300)s. Much 
longer values of r for gh diffusion indicate that it remain 
active for a longer time due to availability of large volume 
of gh. Obtained values of Dr follows Arrhenius behavior 
given by: 

Dr = Doexp{-E/kBT) (4) 

Here, Dq denotes pre-exponential factor, E activa¬ 
tion energy, T temperature and /cp Boltzmann’s con¬ 
stant. The obtained values of E for Dy of Fe and N 
are 1.0±0.2eV and 1.4±0.2eV, respectively while for Fe 
Dgi) it is 0.6±0.2eV. As expected, smaller value of E 
for the gh diffusion signifies that it starts at a relatively 
lower temperature as compared to Dy . Interestingly, we 
observe a slightly higher value of the activation energy 
for Dy of N as compared to Fe (also r for N Dy was 
larger than Fe). Observed discrepancy within Dy of 
Fe and N can be understood due to a stronger metal- 
nitrogen (than metal-metal) interaction predicted theo¬ 
retically for TMMN^’^^. Unlike oxides, TMN are stabi¬ 
lized by a dominant metal-metal interaction and there- 
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fore metal-metal bonds are stronger. However monon- 
itride having ZnS-type or NaCl-type structure, a vol¬ 
ume expansion of fee lattice takes place to accommo¬ 
date N atoms in the interstitial positions. Early theoret¬ 
ical calculations predicted that the interaction distances 
are larger and bond energies are significantly smaller for 
metal-metal bonds as compared to metal-nitrogen bonds 
for 3d TMMN.^^ More recent band structure calculations 
based of density-functional theory, observed localization 
of metal valance bonds leading to weakened metal-metal 
bonding in mononitrides.In additions a suppression 
in N migration vacancy was predicted for slightly off- 
stoichiometry mononitrides of various From 

these theoretical studies it can be inferred that metal- 
metal bonds are readily broken. This situation may lead 
to segregation of Fe atoms in the gh region, leading to 
faster Fe diffusion as compared to N diffusion, under iden¬ 
tical thermodynamic conditions. 


IV. CONCLUSION 

Iron mononitride thin films studied in this work, show 
an anomalous self-diffusion behavior in which N atoms 
tend to diffuse slower than Fe atoms. SIMS depth-profile 


measurements reveal that the diffusion mechanism for Fe 
and N is different. While N diffuses only via a volume- 
type diffusion process, Fe in addition to volume, pre¬ 
dominantly diffuses through grain-boundaries. Fe grain¬ 
boundary diffusion was about two orders of magnitude 
more than its volume diffusion. Even for volume diffu¬ 
sion, N diffusion was found to be less than Fe. This can 
be understood in terms of stronger metal-nitrogen bonds 
(than metal-metal) predicted theoretically for transition 
metal mononitrides, evidenced experimentally in this 
work. 
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